Even-odd effects in magnetoresistance of ferromagnetic domain walls 
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Difference in density of states for the spin's majority and minority bands in a ferromagnet changes 
the electrostatic potential along the domains, introducing the discontinuities of the potential at do- 
main boundaries. The value of discontinuity oscillates with number of domains. Discontinuity 
depends on the positions of domain walls, their motion or collapse of domain walls in applied mag- 
netic field. Large values of magnetoresistance are explained in terms of spin-accumulation. We 
suggest a new type of domain walls in nanowires of itinerant ferromagnets, in which the magne- 
tization vector changes without rotation. Absence of transverse magnetization components allows 
considerable spin accumulation assuming the spin relaxation length, Ls, is large enough. 

PACS numbers: 73.50.-h, 73.61.-r, 75.60.Ch, 75.70.Pa 
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As it was demonstrated first in [g| |2[] , spin accumula- 
tion effects in the presence of a current flowing between 
ferromagnetic and normal metals may result in a consid- 
erable contribution to a contact's resistivity. This phe- 
nomenon is the key element for the GMR devices with 
the CPP (current perpendicular to the plane) geometry. 
The CPP-GMR effect has been studied theoretically for 
the spin- valve systems for a simple triple layer and multi- 
layered structures in the pioneering work j|. The sub- 
sequent experiments (see e.g. j|, [|, |(|) were in excellent 
agreement with the predictions made in j^] concerning 
the dependence of the resistivity on the width of mag- 
netic and non-magnetic components and the role played 
by spin-relaxation mechanisms. Most recently it was dis- 
covered that nanocontacts [Q and domain walls in mag- 
netic nanowires Q possess significant magnetoresistance. 

Realization of different experimental configurations al- 
lows determination of parameters present in the expres- 
sions in Ref. || such as the values of resistivity for each 
of the GMR components and what is most important the 
spin relaxation length, Ls, characterizing the width of a 
non-equilibrium distribution of spins near the contacts 

The formulas in Ref. ||], however, have been derived 
in assumption that while the conductivities of the ma- 
jority and minority spins are different, the corresponding 
densities of states (DOS), g a , remain equal. This as- 
sumption is not realistic. In the presentation below we 
address this issue to demonstrate that taking difference 
in the DOS into account, the changes in the expressions 
for the distribution of the electrostatic potential lead to 
some new observable effects. In Ref. ]?], ||] it was specu- 
lated that the pronounced GMR effects are caused by the 
significant role of spin accumulation. We suggest, as we 



believe, for the first time, that large magnetoresistance 
observed in Ref. || is due to the non-rotational charac- 
ter of the domain walls || which are possible in itinerant 
ferromagnets Eo| . 

Following |3|], we re- write the expression for the cur- 
rent j a = eD a n' a — a a U'{x), (U(x) - the electrostatic 
potential, the index, a = ±, stands for the majority and 
minority spins correspondingly), into the form: 



Ja = — v/x a , 

e 



where 



Ma 



eU{x) 



(1) 



(2) 



is a non-equilibrium (in the presence of a total current, J) 
electrochemical potential for each spin component, and 
the relation D — a/ (e 2 g) is used in ([1], |^). 
The equations for the current components: 

div j a = — (3) 

TS 

together with (0) and the electro-neutrality condition: 



(4) 



present the complete system of equations for each side of 
an interface (in (0) ts is a spin relaxation time). 

To simplify the analysis, we at first assume the ballistic 
regime for the interface, i.e. the width of the correspond- 
ing domain wall is small (contributions due to the spin 
scattering inside the barrier will be discussed at the end 
of the paper). Correspondingly /i Q and j a are taken con- 
tinuous at the boundary. 

In convenient notations of Ref. [0 with the help of 
(0) we obtain for u s = /i + — a_: 
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Subtracting one of Eqs. (|^) from another and making 
use of (H), one obtains the equation for the distribution 
of /is in the "bulk" on each side of the contact: 



V S = —72—: 



where 




(6) 



Equation (^J) coincides in its form with Eq. (14) from 
Ref. H . For our purposes it is convenient to re- write (|^) 



f x"s-» s /L 2 s = S(x)A((x' s ). (6') 

Equation (|^) applies now everywhere in the sample. In 
(^|), S(x) stands for the Dirac delta- function at the 
boundary (x — 0) (equation ((|6^[) is valid when the 
boundary between the domains is abrupt). While us 
is continuous, fi' s has a jump at the interface: 



A(//s) = |u's(+0)-m's(-0)- 



(8) 



The jump can be expressed in terms of the total current, 
J. To do so, we write down the expressions for J, spin 
current js = j+ — j- , to obtain with the help of (|^J^) : 




v+9+ 
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l ,' s (x)-(a + -a.)U'(x). 



(9) 

We specify that on the left from the domain wall spins up 
belong to the majority band. Magnetization changes sign 
across the wall and so does the band's occupation so that 
one has to interchange between <J±,g± 1 in (JsJ) . Since the 
currents are continuous, for the jump of electrochemical 
potential gradient, A(/ig), we get: 



For a multilayer system, Eq. 



J = 
reads: 



C. 



(10) 



n"s - i£ = C]T(-1)*+M(x - an), (i = 1, AT). (11) 
s i=l 

It is also assumed below that the width of the very left 
and very right banks are larger then Ls- The solution 
of ( |TT| ) is a superposition of solutions for a single domain 
wall: 
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(1>S{X) = 



CL S " 



3=1 
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(12) 



We consider first the behavior of the electrostatic po- 
tential, U(x), close to a single domain wall. With the 
help of Eq. (0) one can write: 



g- 



9- 



Hs{x) - 2eU(x). (13) 



Making use of continuity of the electrochemical potential 
at the boundary, one immediately sees from (|l^) that 
potential U(x) is discontinuous with the jump AUi(xi) = 
U[x = Xi, R) — U(x = Xi, L) at the ith domain wall: 



-iy 



9- 



»s(xi). (14) 



The spatial distribution of U (x) can be found by integrat- 
ing the first equation (S) along each side of the domains: 



U'(x) = - 



J 



( g+ 



\9+ 



9- 
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Vis. 



(15) 



The total potential drop across an isolated single do- 
main wall in the chosen geometry is the sum of two terms: 



AU = AUi 



[i s (x)dx 




(16) 



2(T + cr_ (a + + a J) 



J 



(the second term in (|T^) comes about from the currents 
distribution (Q) on the distances of the order of Ls on 
the two sides of the wall) . In the following notations 

<r ± = |(l±/3), . 9± = |(1±<5), 



one obtains: 



AR = 



AU 



= 2 



(17) 



i.e. g± drop out from the total magnetoresistance. Equa- 
tions ( fU3||l7| ) coincide with the results obtained in 0. 
The differences in the DOS for the minority and ma- 
jority spins lead to an appearance of the discontinuities 
( |l4| ) and changes in the dependence of the electrostatic 
potential U(x) (see ([l5])) along the domains. 

To make a numerical estimate for the total potential 
drop we will use the data from Ref. || obtained for Co 
nanowires: Ls — 60nm, p ~ 1.3 • 10 fl-cm. The typical 
values of S are of the order of 0.4 ~ 0.5 ||, |), while f3 is 
of the order of 0.5 -j- 0.7 iQ. After substituting these 
values in (p7|), the resistance drop AR per unit area is 
5 • 10~ n f2, or AR ~ 112 SI for the geometry used in M. 
The value of the potential drop at domain boundary is 



3 
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FIG. 1: Ratio (SR) N /AR is plotted as a function of N @. 
The distance between the domain walls is estimated using the 
formula d ~ \/d w lco (d m ~ lOnm is a width of the domain 
wall, lco — 40nm is the diameter of a wire |8[). 



FIG. 2: Ratio AR co u/AR is plotted as a function of dimen- 
sionless parameter l c — L c /Ls showing the resistivity drop in 
the sample corresponding to the collapse of two domain walls. 
For simplicity, we assumed that there have been four domain 
walls before an external magnetic field was applied. 



approximately the same as the value of the total potential 
drop in ([l6|): 



AXJi/MJ = 8/13. 



(17') 



We would like to emphasize, that the ratio in (17 ) may 
also be negative. 

As an example, let us consider in some more details, 
the drop of the potential, AU(xi), across the very left 
domain for a system of N walls (we also take x\ = 0). 
From (O) and (14) we have 



AU N = 



fa 



a+°-{9+ +9-) 



-JL S S. 



N ■ 



where 
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(-l)V 



-Xj/Ls 



(18) 



(18') 



For the N walls the value of AUn shows an interesting 
"even-odd" effect: 



(-1)" 



(SU) 
fJ 



N = (AU N+1 



AC/, 



N) 



(19) 



Below in Fig. [I] we plot {8R) N /AR = {8U) N /JAR as a 
function of number of the domain walls. The "even-odd" 
effect is well-pronounced at Nd ~ Ls, where d is a size 
of domain. 

Another interesting feature is that changes in Eq. JlS| ) 
could follow a motion of a single domain wall, say at 
x = Xi, through its contribution {—\) i e~ Xi ' Ls into (18'), 



or even a collapse of a domain caused by the applied 
magnetic field (such a collapse has been experimentally 



seen in Ref. |§J). In order to explicitly demonstrate this 
effect, we first introduce the following notations: 



AR C 



N 



AR 



3 = 1 



j/Ls 



(-l) j e- x ^ Ls , 

j=i+l 



The result of our calculations of AR co u / AR as a function 
of the distance between two neighboring domains L c = 
Xi+i — Xi is plotted in Fig. ||. Experimentally, the motion 
of a domain wall can be detected by the STM technique. 

Now we would like to briefly discuss the change in our 
results in case when one takes the finite width of the 
domain wall into account. In particular, we consider the 
importance of depolarizing effects in the Bloch or Neel 
type of domain wall for the results above. Since the width 
of the wall is usually much smaller then the spin diffusion 
length, Lg, the arguments that lead us to results given 
by ( p^jl^ ) still hold. The only modification would come 
from the change in the boundary conditions for the spin 
current. 

Electrons going through the Bloch or Neel domain wall 
loose part of its polarization because the transverse com- 
ponent of magnetization inside the wall creates a torque 
which causes spin's re-orientation. Obviously this process 
will reduce the spin accumulation. If d w is the width of a 
domain wall, spins of electrons traveling through the wall 
are rotated in the transverse component of the exchange 
field, H exc h, in the wall by an angle $ a (a =TI) : 



$ a = (± 
From (prt) one obtains: 



HBHexchd V 



j s (x = +0) - j s (x = -0) = 2en + v 



(20) 



(21) 
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with v given by: 



exch^w 



and after some algebra with help of 
j s (x = +0)-j S (x = -0) = lev- 



9+9- 



9- 



(21') 



■MO)- (22) 



Taking into account (p2|), the expression for the jump of 
li' s at the interface is: 



a(a4) 



j 

<7 + (T- 

2evg + g-{cr + + cr_) 



{a + aJ)(g++gJ) 



(23) 



Ms(0). 



As a result, the expression for the potential drop AU 
across an isolated domain wall acquires the form: 



AU 



J, 



(1-/3 2 ) l + a 
e 2 vgL s <J~ 1 1 - S 2 



(24) 



2 1 - f3 2 ' 
To estimate the value of a, we first re- write the second 



expression in (24) as: 



- 



Ls 
I ' 



(25) 



where I is a mean free path and vf is the Fermi velocity. 
Using the data, provided in |^|, we have: fiBH exc h ~ 
O.leV, d w ~ lOnm, L s /Z ~ 10, u ~ 4.5 ■ 10 7 cm/s, 
«ir ~ 10 8 cm/s. Thus our estimates provide a ~ 14, 
which according to (^3j) significantly reduces the spin 
accumulation effect. This result can also be explained 
slightly differently: from we estimate the depolar- 
ization angle $ ~ 1, i.e. at crossing the Bloch or Neel 
type of domain wall the electronic spins would adiabat- 
ically follow magnetization. To reconcile the above es- 
timates with the significant magnetoresistance experi- 
mentally observed in Ref. ||, |8|, we suggest that the 
domain walls in these experiments were neither of the 
Bloch nor Neel type. Instead, in itinerant ferromagnets 
another type of domain walls, "linear walls", is realized 
|jiof . In a linear wall the direction of magnetization does 
not change, while its absolute value goes through zero in- 
side the wall. Theoretically, "linear" domain walls were 
considered in for local spin ferromagnetic systems at 
temperatures T close to critical temperature Tc U rie- 

To summarize, we have shown that taking into account 
difference in the density of states between the minority 
and majority spin bands drastically changes the distri- 
bution of the electrostatic potential along the domains. 



The discontinuities of the potential across each domain 
wall are of particular interest. The jumps in the values 
of potential can be measured directly by the STM tech- 
nique. Jumps possess such characteristic features as the 
"even-odd" effects counting the total number of domain 
walls in magnetic nanowire in the presence of an exter- 
nal magnetic field. Our results can be directly extended 
to the GMR structures consisting of F/N layers which 
were studied theoretically in || in the approximation of 
equal density of states of majority and minority bands. 
To ascribe large values of magnetoresistance observed in 
|| |8| to spin accumulation effects, it was also necessary to 
suggest that in nanowires made of itinerant ferromagnets 
the domain walls are of a linear type in which magnetiza- 
tion changes without creating a perpendicular component 
that would revert spins of polarized electrons. 
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